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ABSTRACT 

We review Maxwell's multipole vectors, and elucidate some of their mathematical 
properties, with emphasis on the application of this tool to the cosmic microwave 
background (CMB). In particular, for a completely random function on the sphere 
(corresponding to the statistically isotropic Gaussian model of the CMB), we derive 
the full probability density function of the multipole vectors. This function is used to 
analyze the internal configurations of the third-year Wilkinson Microwave Anisotropy 
Probe quadrupole and octopole, and we show the observations are consistent with the 
Gaussian prediction. A particular aspect is the planarity of the octopole, which we 
find not to be anomalous. 
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1 INTRODUCTION 

The cosmic microwave background (CMB) is a unique cos- 
mological probe. It is a window onto the early Universe, a 
backlight on the more recent, and it defines our observable 
Universe. It is therefore our best tool for investigating pre- 
dictions of inflation theories, and cosmological scale proper- 
ties of space, such as statistical isotropy. Predictions from 
inflation, such as homogeneity, isotropy, and the existence 
of Gaussian fluctuations, currently form the basis for our 
fundamental framework, with consequences for all other cos- 
mological studies. It is therefore essential that we test these 
predictions, with appropriate tools. Current and future ob- 
servations of the cosmic microwave background (CMB) pro- 
vide us with a wealth of cosmological information, but they 
also confront us with many challenges. The large data-sets 
require novel data analysis and map-making techniques if 
we are to thoroughly exploit the data of its scientific infor- 
mation. 

Since its launch in 2001, NASA's Wilkinson Microwave 
Anisotropy Probe (WMAP) has observed the CMB over the 
full-sky, in five frequency bands between 22 and 90 GHz. 
Multiple frequency coverage has enabled reliable separation 
of the Galactic foreground signa l from the CMB ani s otropy , 
producing spectacular results (jSpergel et al.l |2003| , l2006h . 
The observations provide further support for the ACDM 
concordance cosmological model, and are remarkably con- 
sistent with a range of independent observations such as 
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the clustering of large-scale structure, and luminosity dis- 
tances from Supernovae observations. However, some inter- 
esting features have been reported, on large-scales, which 
are anomalous with respect to our standard cosmological 
framework of a Gaussian statistically isotropic CMB. If an 
alternative framework is required, then this has far reach- 
ing consequences for the interpretation of all cosmological 
observations, and fo r the determination o f cosmological pa- 
rameter constraints ([Hansen et al ] |2004al ). 

On multipole scales ^ < 60 there is an asymmet- 
rical distribution of power on the CMB sk y, with more 
power in th e southern ecliptic hemisphere (Eriksen et al.l 
|2007l, l2004al : [Hansen et all l2004bl : iDonoghue fc Donoghud 
20051 ). Similar asymmetrical signals have been seen with 
the bispectrum, n-point correlation functions, spheri- 
cal wavelets, Mi nkowski functional s , genus statistic, and 
local curvature (Land fc Magueiid 2005bl: Eriksen et al.l 



2004 lEriksen et al.H2004bl : iParkI l20o3 



20051: IViel va et al 
Han sen et al. 2004^) 

A separate 'anomaly' involves just the largest scales, 
with I = 2,3 co rrelated with eac h other a nd with the 
dipole direction ([Copi et al.l [200i . [2006; Ra lston fc Jaiiil 
l2004l : Ide Oliveira-Costa et al.l |2004| : [Copi et al.l |2007| ) . Fur- 
thermore, I = 4, 5 are also seen to align with this Axis 
of Evil' direction (|Land fc Magueiioll2005ah . Independently, 
^ = 3 are reported to be particularly planar |Tegmark et al.l 
I2OO3I : Ide Oliveira-Costa et all |2004| : ISchwarz et all 120041') 
and the power on these scales is marginally low ( Efstathioul 
[20041 . The multipole vectors, described below, are an im- 
portant tool in these studies. 

However, the interpretation of these results is unclear. 
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Different methods are employed, and various statistical mea- 
sures invented, without systematically examining underly- 
ing mathematical properties of the representation. It is also 
not clear if some 'anomalies' are connected and thus over- 
reported, and a more systematic basis for tests of non- 
Gaussianity of the CMB is required. The usual harmonic 
n-point correlation functions, the bispectrum, trispectrum, 
etc., are unattractive as they are computationally intensive, 
abstract with relation to the underlying properties of the 
map, mutually dependent, contain redundant information, 
and are only meaningful under the assumption of statistical 
isotropy. 

In this paper we revisit the multipole vectors, and study 
their mutual probability distribution for Gaussian maps. 
As an example of the general results, we consider the pla- 
narity of the octopole, but not the possibility of correla- 
tions in direction between different multipoles. The multi- 
pole vectors are a geometric representation of the spheri- 
cal CMB map which is explicitly rotationally invariant, and 
whose probability distributions can be calculated for vari- 
ous Gaussian and non-Gaussian distributions of the CMB. 
They were introduced into cosmology and numerically an- 
alyzed by Copi et al. (|2004|'). and we re subs e quent l y stud - 
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led bv iLand fc Magueiiol (I2005dll3): iDennisI (12004 
iLachieze-RevI l|2004l ): iKatz fc Weeks! (|2004l ) : I Weeki 
among others. More recent astrophysical applications of 
the m ultipole vectors include further analysis of the CMB 
data (jRakic fc Schwar3 |2007[), characterizing the pe culiar 
velocity field of supernovae l|Haugboelle et aU |2006| ) , and 
quantify ing systematic effects of dipole straylight contam- 
ination dCruppuso ct al. 2007*). 

Most of the multipole vector CMB analysis has in- 
volved a comparison of various derived quantities with the 
equivalent values from numerical Monte Carlo simulations; 
here, we present an alternative approach, based on ana- 
lytic evaluation of the multipole vector probability distri- 
butions for statistically isotropic Gaussian random maps. 
We compute the full probability density function for the to- 
tal multipole vector configuration, and use this to predict 
analytically the statistics of the multipole vector configura- 
tion. In Section [2] we introduce the multipole vectors and 
describe their numerical evaluation using polynomial root- 
finding algorithms (through the Majorana polynomial for- 
malism) ; in Section[3]we review their properties for Gaussian 
random maps, and present the analytic Gaussian multipole 
vector probability density (calculated in the Appendix). In 
Section |4l these general statistical results are compared to 
the recent WMAP data at the largest scales, namely the 
quadrupole and octopole. In particular, we will describe an 
analytically tractable measure of planarity for the octopole, 
and demonstrate the WMAP data is very close to the Gaus- 
sian prediction. We conclude with a discussion in Section [5] 
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where m = —£, —£ + 1, —1, 0, 1, £—1,1, and is an 
associated Legendre polynomial. Yem{d,4') = Yim{n) is the 
angular portion of the general solution to Laplace's equation 
in spherical polar coordinates, and £ corresponds to the total 
angular momentum quantum number, and corresponds to 
an angular size of ^ ~ 180° /^. All information in the map 
is now contained in the spherical harmonic coefficients aim, 
with reality of / enforced by the condition 



at.- 



= (-1) 



(3) 



The expansion ([1} for a fixed £ ('multipole') has 2£ + 1 real 
degrees of freedom. 

Following the Copernican principle - that there should 
be no special observers in the Universe - statistical isotropy 
(SI) is a natural assumption to make about the Universe. As- 
suming SI for all observers further implies homogeneity, and 
thus it is common for us to adopt the most general homo- 
geneous and isotropic framework in cosmology, namely the 
Friedmann-Robertson- Walker metric. However, alternative 
frameworks exist, such as the anisotropic Bianchi models 
or the inhomogeneous Lemaitre-Tolman-Bondi models. Pre- 
ferred directions would also exist if the Universe had a non- 
trivial topology, or contained large-scale magnetic fields. SI 
of the CMB map is equivalent to the property that no quan- 
tity depends on any particular frame of reference - the map 
is statistically invariant to global rotations. Deviation from 
SI in the CMB observations would imply exotic physics in 
the early Universe or, perhaps more realistically, residual 
contamination by anisotropic foregrounds (e.g. emissions by 
Galactic diffuse components) or by residual systematic ef- 
fects. 

Such statistical properties of the spherical map can eas- 
ily be conferred to the a^mS. We assume that any cosmic 
probability distribution of the a^^s is ergodic, that is, spa- 
tial averaging is equivalent to ensemble averaging, and such 
averaging is denoted {•). Therefore, the 2-point correlation 
function can be written in terms of the aim coefficients, 

(/("^/("•a)) = ^ {aiimiae2m2)Yi^miini)Ye2m2{n2), 

^1,2, "11, 2 

(4) 

and the action of a rotation R on the spherical harmonics 
involves the Wigner O-function (rotation matrix element) 



(5) 



The orthogonality condition of the ©-function, 

y ^ '^rnim.'D-m2-m ~ Smim, (6) 



2 THE MULTIPOLE VECTORS 

It is natural to represent data on the direction sphere (such 
as the cosmic microwave background) through the Fourier 
expansion 



(1) 



implies that {/(ni)/(n2)) is rotationally invariant (i.e. sta- 
tistically isotropic) if and only if 



(a^jmidfama) — 3l-^l2Sm.im2Cei, 



(7) 



where Ci is the angular power spectrum. 

The issue of Gaussianity is separate from SI, and if the 
CMB fluctuations are Gaussian, then all the information in 
the map is contained in the 2-point correlation function, or 
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equivalently in the angular power spectrum Ce as defined 
in ((TJ; higher order functions - which otherwise probe non- 
Gaussianity - depend only on the 2-point correlation func- 
tion, or are zero. SI is a consequence of Gaussianity, but 
there are other, non-Gaussian, statistically isotropic ensem- 
bles; this will be discussed further in Section [3] 

An alternative representation to of a function on a 
sphere is via the multipole vectors (MVs). For a real multi- 
pole of order I, represented 

t 

Me,(t>)^ ae^Yi^{e,(l}), (8) 

with coefficients satisfying ([Sjl (i.e. fe{x,y,z) is an eigen- 
function of the Laplacian on the unit sphere [x^ + + = 
1) with eigenvalue —£{£ + 1)), there exist £ unit vectors 
vi, . . . ,ve such that 

feix,y,z)^AD^,...D^,^, (9) 

where D^. = • V is a directional derivative operator, r = 
{x^ + + z'^y^^ and j4 is a real co nstant. A mor e useful 
form of the representation is given bv lDennid (|2004| ). 

fe{n) = B{n-vi)...{n-ve)+r^T, (10) 

where B is another real constant, and T is a scalar, given by 
the contraction of a certain Cartesian tensor of rank £ — 2 
with unit vector n, which is fully determined by the 2£ com- 
ponents of the Vi and is completely traceless and symmet- 
ric (Zou & Zheng 2003). A graphical representation of the 
^ = 2, 3, 4 components of the WMAP temperature data of 
the CMB, together with their multipole vectors, is shown in 
Figure [T] 

A simple way of connecting the two multipole repre- 
sentations (|8]) and (|10p is using stereographic projection 
of the direction sphere (as the Riemann sphere) into the 
complex plane; thus the point on the sphere with coordi- 
nates 9, (j) is stereographically mapped to the complex num- 
ber e*'* tan(6'/2). The stereographic images of the multipole 
vectors appear as th e roots of th e Majoran a polynomial 

i spin coherent state) (|Dennisll2004l ) (also see iHelling et al.] 
200i)), 

t. 

P{0 = ^ a.lmlJ-lmC"^\ (11) 

where ^iim is a numerical factor. 

Thus p{C,) is the complex polynomial acquired by replacing, 
in the expression ((Sjl, Y^m with /J^mC^^™- By the fundamen- 
tal theorem of algebra, piC,) can be factorized: 

21 I 

PiO = - = - ^»)(^ + VC*)- (13) 

■i=l i-1 

The vectors corresponding to the roots Q form antipodal 
pairs since ae-m = (— l)'"a£„, implying that if C,i is a root, 
so is its antipode on the Riemann sphere — 1/C* (jPennisI 
|2004 ). and we label the roots such that (^i^e = —l/Q, i.e. 
(iXi+e are antipodal. The roots can be written in terms of 
modulus and argument, 

=tan(6li/2)exp(i<^i)> (14) 




Figure 1. Representation of the multipole vectors for the 
WMAP3 multipoles £ = 2,3, 4, 5. The colour map on the spheres 
represents the underlying WMAP temperature data. 

with 9i , the spherical coordinates of the ith multipole 
vector. 

The representations (|11[) . (|13|l provide a simple method 
for computing the MVs, since the problem is reduced to 
finding the complex roots of a polyn omial of or d er 2£ . 
Alternative methods include those of ICopi et al.l (|2004D 
which use an iterative algorithm to solve (|10p . and a simi- 
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lar ap proach was previously outlined bv lMezei fc Campbelll 
l| 19761 ) ■ albeit in a different physical context. Obviously, 
all methods have limited numerical accuracy, due to the 
global nature of the problem. We note that numerical root- 
finding algorithms run into fundamental difficulties for high- 
order polynomials, due to problems of ill-conditioning and 
the limited numerical accuracy of the polynomial coeffi- 
cients (i.e. the aim coefficients); with 2 decimal place ac- 
curacy in the coefficient s, practical r oot-finding is proba- 
bly limited to ^ < 20 (| ActonI [l970l ) . It is harder to lo- 
cate the source of numerical error in the iterative ap- 
proach. On account of these numerical problems, MVs have 
been studied in the CMB for lower multipoles, particu- 



p^ 

larly for £ ^ 10 (ICopi et al.ll2004l: iLand fc Magueiioll2005d 
Weeks! l2004l: iRakic fc Schwarj|2007l: ICruppuso et all 12007 



Haugboelle et al.ll2006l : ICopi et al.ll2006l ) 

The MVs are, in fact, directors ('headless vectors'), 
since their sign is globally undefined (a ± transform on any 
pair of Vi leaves ((9)1, (|10|l unchanged). Furthermore, just like 
roots of a polynomial, they are not ordered amongst them- 
selves. Counting the degrees of freedom we see that for a 
multipole £, there are 2£ degrees of freedom in the £ unit 
vectors, and one in the constant ^4 in © (or B in (llOp . 
agreeing with the 2^-1-1 freedoms in the airn coefficients) . 

A simple example of the multipole vectors is those of a 
pure, real spherical harmonic KeYim- Of the £ MVs, £— \m\ 
are aligned with the z-axis (in the Yim frame) and the re- 
maining \m\ vectors lie in the a;j/-plane, with vertices forming 
a regular 2m-polygon. Since this system of vectors rotates 
with the function, it is easy to identify pure harmonic modes 
in any frame of reference; the distribution of the a^mS of a 
pure mode in an arbitrary frame of reference is certainly not 
so transparent. This also applies to arbitrary multipoles, 
whose configuration of MVs rotates directly with rotation 
of the function, unlike the coefficients aim (which change 
according to the appropriate ©-functions). This property 
makes the MVs an extremely useful tool for studying the sta- 
tistical isotropy of the CMB - their configuration and statis- 
tics do not depend on the frame of reference. This avoids the 
problems inherent in manipulating the a^mS directly, and of 
having to find complicated rotationally invariant functions 
(bispectrum, trispectrum, etc.) which are hard to connect 
to underlying properties of the map and increasingly com- 
putationally intensive to compute. A drawback of the MVs 
is the problem of their numerical determination, described 
above. 

The multipole vector formalism was apparently first 
introduced by Maxwell (jMaxwcU 1 89ll ), an d is discussed 
in several classic mathe matical textbooks llHobson) Il93ll: 
ICourant fc Hilbertlll953l ) (also see the references in lOennisI 
It has been appl i ed in studies of quantum chem- 
istry (^Mozoi fc Campbell' 'l97^) and continuum mechan- 
ics (Backus 1970; Zou fc Zheng 2003). Its generalization to 
arbitrary spin states (where £ may be a half-integer and 
the coefficients a^™ do not necessarily satisfy Q) leads to 
the imp ortant canonical decompo sition of totally symmetric 
spinors (jPenrose fc Rindleiill984 ) and the stellar representa- 
tion ofs£in_coherent_^ (Bacry 
1 19741 : iBogomolnv et all 1 19961 : ISchupp|fl999f) . In cosmology, 
the MV method has received considerable attentio n, from 
mathematicians exploring features of the formalism (iDennid 
I2OO5I : iKatz fc Wee"k3l2004l : lLachieze-Revll2004l : IWeeksll2004 



iHelfing et al.l I2OO6I ) , and from ast rophysicists applying it 
to the WMAP data of the CMB l|Copi et all |2006| . 12004 
iLand fc Magueiidl2005dl ). 

In the following section, we will discuss the implications 
of different statistical assumptions of the aemS on the prob- 
ability distribution of the multipole vectors. 



3 GAUSSIAN MODEL 

A candidate ensemble for the atm coefficients is 'Gaussian', 
that is, they are independent, identically distributed com- 
plex Gaussian random variables. Gaussianity of the CMB 
fluctuations is a generic prediction of simple inflation the- 
ories, because during slow-roll inflation the near flatness of 
the potential results in each Fourier mode of the inflaton 
perturbation evolving independently. This property trans- 
lates to the density perturbations during reheating, and the 
central limit theorem implies their sum (i.e. into real space) 
will be Gaussian. Furthermore, the probability distribution 
for each Fourier mode is already expected to be Gaussian (if 
the perturbation is in the vacuum state). These models also 
predict that the Fourier modes have a near scale-invariant 
power spectrum. Even in these simple models, deviations 
from Gaussianity are expected on small scales where the 
central limit t heorem breaks down. More complica ted in- 
flation models l|Wandsll2007l : iRigopoulos e"rZll2006l ). or al- 
ternative t heories of stru cture formation such as topologi- 
cal defects (|Battvelll999l ). can predict higher levels of non- 
Gaussianity. 

Requiring SI alone is a far weaker assumption, that 
the probability density function (PDF) for the set of a^mS 
depends on rotationally invariant quantities (the angular 
power spectrum Ce, bispectrum, trispectrum, etc.), and does 
not fix the ensemble uniquely. As described in the previous 
section, any probability density function (PDF) for the set 
of aim coefficients satisfying ((T} is statistically isotropic. For 
instance, it is possible to construct statistically isotropic en- 
sembles which have a fixed MV configuration that is ran- 
domly oriented (for instance, for £ — 2, the two MVs have 
the same - random - direction, equivalent to Y20 in some 
random frame). A Gaussian PDF, in which the atrnS are in- 
dependent and Gaussian (with identical variance for same 
£) clearly satisfies ((T}, and is unique in the property that 
its 2-point correlation function (or equivalently, its power 
spectrum) determines the ensemble uniquely; statistically 
isotropic non- Gaussian ensembles, whilst satisfying have 
different n-point correlation functions (equivalently, aver- 
ages of more than two airn coefficients are non-Gaussian). 

Nevertheless, there are other statistically isotropic 
PDFs which have the same MV distribution as the Gaus- 
sian model - using the MVs alone, it is impossible to dis- 
tinguish Gaussianity uniquely from certain other ensem- 
bles. This wider class of ensembles is characterized by the 
property that the only quantity appearing in the PDF of 



the aimS is Ce 



\ae„ 



/{2£ + 1); equivalently, the 



{2£ + l)-dimensional vector {ai^-i, an) is invariant to 
all {2£ + l)-dimensional unitary transformations preserv- 
ing ([Sjl. In iDenni^ (|2005l ). this wider class of ensembles was 
called completely random, as the statistics of the multipole 
ft depend only on the distribution of its normalization Ce, 
and the MVs have the same statistics for any completely 
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random ensemble, independent of the precise PDF of Ci. 
This may be understood from the polynomial representa- 
tion the MVs are distributed by the polynomial roots, 
whose distribution is independent of the any overall mul- 
tiplicative factor of p(C)) equivalent to Ci. Clearly, Gaus- 
sian distribution of the atmS satisfies the condition to be 
completely random, and the PDF of Ci is that of the ra- 
dius in a multidimensional isotropic Laplace distribution. 
A non- Gaussian example of a completely random ensem- 
ble is that in which Ci is absolutely fixed (its PDF is a 
(5-function); of course, such a distribution does not have 
any cosmological significance. Another completely random 
example is that of Gaussian and independent a^^s, apart 
from a cut forcing a maximum Ci. Such an ensemble was 
explored bv lRakic fc Schward (|2007h - indeed, they observe 
that in this case the MVs are indistinguishable from the 
independent Gaussian case. If, in addition to complete ran- 
domness, the aimS. are also assumed independent (i.e. the 
full PDF factorizes, P{{airn}) = P{ae,-i)---P{ae^i)), then 
the completely random ensemble must be Gaussian (the ar- 
gument is analogous to that of the Maxwell distribution in 
kinetic theory (Dennis 2005)). Another feature of all com- 
pletely random ensembles - including the Gaussian - is that 
the PDFs for separate multipoles are completely indepen- 
dent. We will not examine this prediction against the CMB 
data in the present paper. 

Perhaps counterintuitively, for a given i, the MVs in the 
completely random case are not statistically independent, 
but have a unique joint probability distribution, which we 
derive exactly below. The derivation uses the fact that the 
MV distribution is equivalent to the probability distribution 
of the roots of the random Majorana polynomial and is 
a special case of more gener a l investigations of random poly- 
nomi als (e.g. lHannavlll996l : iBogomolnv et ahl Il996l : iProsenI 
1 19961 ). where independent Gaussian random coefficients are 
assumed, although not the condition (|3]). In our analysis in 
Section[4l we will compare the MV distributions of the com- 
pletely random model (including the Gaussian model) with 
that of independent MVs; although this latter does not have 
any obvious cosmological application, it is a non- Gaussian 
ensemble satisfying SI whose statistics are very simple to 
calculate analytically. 

It has been pre viously shown ()Dennisl l2005l : 
iLand fc Magueiid l2005d ) that the 2-point function of 
the MVs exhibits quadratic repulsion (as usual for the roots 
of ra ndom polynomials and matrices (jBoeomolnv et aD 
1 19961 )). Given this repulsion, one may expect that the most 
probable MV distribution for fixed £ is at the vertices of the 
regular polyhedron with 2£ vertices, if it exists. If it does 
not exist, then there will be no unique maximal probability 
distribution (similar constructions fo r polynomials on 
the R iemann sphere are considered by lAtivah fc Sutcliffa 

For any completely random distribution of the aim co- 
efficients, including the Gaussian distribution, the probabil- 
ity distribution for the £ multipole vectors vi, ...,vi can be 
found analytically, based on a transformation of the set of 
Gaussian-distributed coefficients {aim} to the set of £ inde- 
pendent roots {(i} of the Majorana polynomial (lll|) . The full 
details of this calculation are given in the Appendix. For the 
multipole of order £, the full joint probability density func- 
tion of the £ multipole vectors {vi}, assuming dependence 



only on Ci, is given by 
{2£ 



l)!!n-iij! 



(2ny2^'£\ 

n 



(15) 



3«a{{.a,<.)/2 Ylf^^ cos 



e+1/2 ' 



where Ojk denotes the angle between the MVs labelled by j 
and fc, and in the denominator, summation is over permu- 
tations a in the set S2i of permutations on the 2£ indices 
j = 1, ...,2^; such a permutation gives rise to a system of 
piecewise-geodesic curves on the direction sphere (as each 
Vj goes to ^^{j) to t'a-2(j), ••■)) ^-nd is the signed 

solid angle on the sphere enclosed by these curves. It should 
be noted that vi+i = —Vi, and the term in the denominator 
involves the 2£-element set {vi, ...,V2e}- 

In practice, it is easier to work directly with the roots of 
the Majorana polynomial (stereographic images of the MVs 
in the complex plane), which have the probability distribu- 
tion equivalent to (fTS)) . 



{2£- 



l)!!n-lij! 



{27vy£iuU\CA' 



n 



2e 

j,k — l,j<k 



e+i/2 ■ 



(16) 



The function p5|) . and especially its stereographic coun- 
terpart (|16p . bears some similarity to the full probability 
distributions of roots of ot her Gaussian random polynomi- 
als re presenting spin states (IBogomolnv et al ] |l996l : lHannavl 
Il996l ): the discriminant of the polynomial in the numerator 
is divided by a combinatorial factor which mixes the roots 
nontrivially. 

The functions Pi{{vi}) and Pe{{C,i})) contain all the 
statistical information of the directions of the MVs, and the 
functions for different £ are mutually independent; the func- 
tion (|15|) is one of the main results of this paper. The ex- 
pressions (|15p and (|16p are unwieldy for high £, primarily 
because the sum in the denominator involves (2£)\ terms, 
most o f which are no ntrivial. 

In lDenn5 (12005V the 2-point multipole vector functions 
pi were calculated (corresponding to Pi{{C,i}) only when £ — 
2). Normalized with respect to angle 9 (between and 7r/2), 
this is 



pi{e) = sine (^(£- 1)(1 + cos6')^_D 



5/2 



{{2£D - 4buv - (6^ + v^){a - 1 - u^)) x 
{dD — 2cuv{a -j- 1 — u^) — {c' + av'^){a — 1 — u^)) 
+ {2£D - 2cuv - buv{a + 1 - m^) - u^(a - 1 + m^) 
—bc{a — 1 — it^))^ + {wD — 2bcu — uv'^{a + 1 — u^) 
-bv{a-l + u'^)'-cv{a~l-u'^)f} , (17) 

with D — {a — 1 — — 2m) (a — 1 ~ + 2u) and 



a 



sec 

d = 2^sec«-2 e 



2, b — 2£tan j , c = , 



sm o sec 



M 



2 ^1 -f 2£tanf) 



-2/tan 



21-1 



-2£(2£ 



u = tan 
- 1) tan 



2 ' 
21 9 
2 ' 

2^-2 e 



(18) 



For the remainder of the paper, we will compare the full 
joint probability Pi and the 2-point pe for the lowest two 
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nontrivial ^-components of the CMB distribution, namely 
the quadrupole £ — 2, and octopole £ = 3. This probability 
is shared by the Gaussian and other completely random en- 
sembles; we will also compare the results for the ensemble 
of independent MVs. We will not be concerned here with is- 
sues of data analysis and error in determining the multipole 
vectors, but merely note that the full PDF enables the er- 
ror in the multipole vector directions to be mapped directly 
to t hose in the aemS- The £ = 2 case wa s considered be- 
fore l|Land fc Magueiidl2005j : lDennisll2005l ). but is included 
here for completeness; the £ — 3 analysis is new. 



4 COMPARISON WITH CMB DATA, 



2,3 



For tests of statistical isotropy it is essential that we use full- 
sky maps, as a mask induces a relatively large anisotropic 
signal of its own. We therefore examine the 'internal-linear- 
combination' (I LC) maps of the WM AP collaboration. The 
WMAP mission l|Spergel et al.ll2003l ) produced full sky CMB 
maps from ten differencing assemblies (DAs). They also pro- 
duced an ILC map, from combining smoothed frequency 
maps with weights chosen to minimize the rms fluctuations, 
using separate sets of weights for 12 disjoint sky regions. 
In the first-year data release the WMAP collaboration ad- 
vised that the ILC map be used only as a visual tool. How- 
ever, for the third-year release a thorough error analysis of 
the ILC map was performed, and a bias correction imple- 
mented ( Hirishaw et al.. 20061 ). The resulting third-year ILC 
map (herein WMAP3) is expected to be clean enough on 
scales ^ < 10 to be used without a maslfl. 

As a consistency check we compare the WMAP3 data 
with some third-party maps, in which different methods of 
removing the for eground signal were e mployed. The third- 
year ILC map of iTegmark et al.l l|2003l ) (TOH herein), cal- 
culated using an 'internal' method (as WMAP), but with 
weights depending on scale in harmonic space as well 
as Galactic latitude. This is advantageous since different 
sources of contamination dominate at different scales - fore- 
grounds at large scales , and noise at smaller sc ales. We 
also examine the map of jpark. Park, fc GottI (|2007l ') (PARK 
herein), who employed an ILC method but fitting differ- 
ent weights over hundreds of disjoint regi ons rather than 
j ust 1 2. The final datset is that produced bv lBielewicz et al] 
l|2005h . who reconstructed the low-multipoles from high- 
latitude WMAP data by means of the power equilization 
filter method. We use their V-band map (PE herein). 

For £ — 2, in the Galactic frame, the multipole vectors 
of the WMAP3 are 



/ -0.562 \ 


/ 0.971 \ 




0.815 


0.048 


(19) 


\ 0.138 / 


\ 0.234 / 





The angle between the quadrupole CMB vectors is 61.7° 
(to 1 decimal place). The resulting vectors are very simi- 
lar for our third-party maps, although the inter-angle varies 
somewhat; {TOH, PARK, PE} = {82.8°, 74.5°, 53.1°}. Due 
to the low number of degrees of freedom, the low-^ multi- 
poles are particul arly sensitive to the various tr eatments of 
the galactic plane (jNaselskv. Verkhodanov. fc Nielsen 2007) . 

^ WMAP data available from Ihttp: / /lambda. gsfc.nasa. gov] 



However, as we will see, this variation does not effect our 
conclusions. Choosing the MVs (rather than their antipodes) 
such that the z-component of the MV is positive, the £ — 3 
WMAP3 vectors are 



0.905 
0.401 
0.145 



0.677 
-0.708 
0.202 



-0.051 
0.770 
0.636 



(20) 



We label these ito,wi,W2 in sequence, and write down the 
three cosines between these vectors (which define the con- 
figuration uniquely) 

uo ■ ui — 0.358, ui ■ U2 = —0.451, W2 • uo = 0.353, (21) 

corresponding to angles between the directors of 
69.0°, 63.2°, 69.3° respectively. The MV director angles for 
the third-party maps are 



TOH 
PARK 
PE 



68.9°, 66.9°, 73.8° 
65.4°, 64.4°, 73.5° 
64.5°, 62.1°, 65.4° 



(22) 



We now analyze these MVs with the Gaussian predicted 
statistics of Section [3] In subsection 14.11 we will discuss the 
2-point MV correlation function and its comparison with 
the Gaussian 2-point function (I17|) : in subsection 14.21 we 
will discuss how the octopole data compares with P-i{{vi}), 
and in subsection 14.31 we will discuss a natural measure 
of non-planarity for the octopole whose distribution can be 
almost-analytically handled by the Gaussian MV distribu- 
tion function P3. 



4.1 2-point multipole vector correlation functions 

The 2-point correl ation function probabilities were previ- 
o usly calculated in Demii^ (|2005l ). and studied numerically 
in lLand fc Magueiid (|2005cl ). When £^2, this is 



P2(cos6') = 



27(1 - cos^ 61) 



(23) 



(3 + cos2 6')V2 ■ 

Since the angle 8 between the two multipole axes is the 
only degree of freedom for the £ = 2 configuration, this is 
also the joint PDF P^dvijV^}) of the multipole vectors 
for the quadrupole, and can be explicitly shown from (|15|) . 
Normalizing over 0^9;^ 7r/2, therefore, we have 



P2i0) 



27sin'^6' 



(24) 



(3 + COs2 61)5/2 • 

In Fig. [2ji, P2 is plotted together with the WMAP3 
data SwMAPS of 61.7°. The most preferred angle is 90°, with 
a probability of Prob(S > 80°) = 0.29, and a median of 
9 « 72.2°. For the WMAP3 data we find Prob(6l < 61.7°) = 
0.28, giving no indication of any significant discord between 
the data and the Gaussian (completely random) prediction. 
Similarly, the results from the third-party maps give are 
not significant as the lowest angle (for PE) finds Prob(5 < 
53.1°) = 0.16. Thus there appears to be nothing anomalous 
about the distribution of the multipole vectors for £ = 2. 

When £ — 3, the 2-point function is, from (|17p . 



P3(A) = 



75^3(1 - X^)(3 + A^)^(88 + 67X^ + 34X^ + 3X'^ 
{{■i-X + A2)(4 + X + X2)(5 + 3A2))5/2 

(25) 

with X — cos 9. As with P2, the maximum occurs when 
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Figure 2. The 2-point correlation functions for Gaussian ran- 
dom multipoles, and the CMB WMAP3 data, a) £ = 2, with 
angular separations 6wMAP3 = 61.7°; b) ^ = 3, with 6wmap3 = 
63.2°, 69.0°, 69.3°. 



9 — 90° - when the vectors are orthogonal. The function is 
plotted in Fig.[2]D, together with the data from (|2ip . We find, 
for the CMB data values, Prob(6' < 6*3) = 0.29, 0.41, 0.42 for 
63 = 63.2, 69.0, 69.3° respectively, with very similar results 
for the third-party maps. Thus, again the 2-point MV corre- 
lation function does not highlight anything anomalous about 
the octopole, although we look to the full joint PDF consid- 
ered in the next subsection for a more complete analysis. 



4.2 Full octopole Gaussian PDF 

Since the PDF of any statistically isotropic ensemble is in- 
variant to the choice of reference frame for the MVs, it is 
convenient to rotate the system of MVs vo, fi, f 2 so that vo 
points in the -j-z-direction and vi lies in the x, 2-plane 

There are therefore three parameters which govern the 
internal geometry of the system: the z-coordinate ci oi vi 
(i.e. Ci — cos^i), the 2-coordinate C2 of V2 (i.e. C2 — cos ^2) 
and the relative azimuthal angle (p between vi and V2- Our 
chosen standard coordinate system is therefore 




Sl 



Cl 



V2 




(26) 



where Sj — sinBj = 



1 — c~,j = 1, 2. However, due to the 
antipodal nature of the MVs, it is not appropriate to define 
Cl , C2 over the whole sphere, and according to the labelling 
convention of (|20p . ci, C2 J? (i.e. the MV axis is labelled by 
the vector in the northern hemisphere). 

Using (|21[) . it is straightforward to determine that, for 
the CMB MV data (with the same ordering), 

Cl,WMAP3 = Uo ■ Ui = 0.358, 

C2,WMAP3 = Uo • 112 = 0.353, (27) 
(/>WMAP3 = arg(uo X ui) ■ {uo x U2 + 2112) = -2.293. 

These are readily checked by verifying that the scalar prod- 
ucts of (|26p agree with (|21[) . and that the scalar triple prod- 
uct is the same for the two sets of vectors. 

The full £ = 3 octopole PDF P3 is given by ^ (ac- 
tually, computed using (|16|l ). is significantly simplified by 
choosing the MVs in the standard coordinate system (I26|) . 



reducing to 

-P3(C1,C2,0) = 

750V3s?si(l 



(28) 



(ciC2 + S1S2 COS ( 



7r(5 + 3c^ + 3c| + (ciC2 + S1S2 cos (j)){ciC2 + 3siS2 cos 4>)y^^ 

It is evident that Ci and C2 appear on an equal footing, 
as they should. Further, because the multipole vectors are 
not naturally ordered, this PDF does not depend on their 
ordering (although the exact equation for parameter (p in 
H27|) depends on the xy-quadrant U2 occupies in the rotated 
standard coordinate system). 
P3 is normalized, i.e. 



1 .1 
dci / ciC2 
Jo 



)P3(ci,C2,0) = 1, 



(29) 



which is easily checked numerically. Our choice of parame- 
ters is such that the probability density element dci dc2 dij) 
is uniform. In the model in which the MVs are statistically 
isotropic and completely independent, the PDF for ci,C2,<^ 
(over the same integration domain) would evidently be con- 
stant l/27r. 

The 3-dimensional distribution of P3(ci, C2, (j)) is shown 
in Fig. |3] It is clear from the figure that the maximally prob- 
ably distribution of the MV axes is mutually perpendicular, 
that is Cl = C2 = 0, = ±7r/2, or \vo ■ (vi x V2)\ = 1. 
This supports the previous statement that the unique max- 
imal probability configuration exists when 2£ is the num- 
ber of vertices of a regular polyhedron, in this case an oc- 
tahedron (whose six vertices are mutually perpendicular), 
and the MV configuration falls into the octahedral symme- 
try group (the octopole function with maximum probability 
(oc xyz oc iYs — iY^^) has tetrahedral symmetry, not being 
inversion-symmetric) . 

This maximal probability occurs with probability den- 
sity 6\/3/\/57f ~ 1.48, about 9 times that for independent 
MVs. The probability is low when the MVs are nearby, 
demonstrating repulsion as with the 2-point MV function. 

The MV directions of the octopole CMB data are rep- 
resented in Fig. 31 and this has probability density 0.362, 
which can be compared to the maximum (orthogonal con- 
figuration) with a probability density of 1.479, and the min- 
imum (all aligned) with zero. We find that 87.4% of the 
{ci, C2, 0} configuration space has a probability density less 
than 0.362, and thus the WMAP3 octopole MV configura- 
tion is not improbable. The third-party maps (TOH, PARK, 
PE) have octopole configurations with probability densities 
0.484, 0.380, 0.252 respectively (with 91.5%, 88.2%, 81.2% of 
configuration space with lower probability density). Thus, 
according to the Gaussian multipole vector joint probabil- 
ity density function, there is nothing particularly anomalous 
about the £ = 3 configuration of the multipole vectors, in 
agreement with iLand fc Magueiid (l2005al . 120061). 



This contrasts with the fi nd ings of iTegmark et al 
(200:|; Ide Oliveira-Costa et al.l (|2004D : ISchwarz et~ 
(2(X)J), in which the octopole data is interpreted to be 
remarkably planar. That is to say, the multipole vectors 
tend to lie equally spaced in a plane, or equivalently the 
octopole is dominated by 13.3 in s ome fr ame. These claims 
were re-examined in ICopi et al.l l|2006l ). where planarity 
was found with a reduced evidence of only 89% Confidence 
Limit, using the same statistic. Indeed the vectors do 
tend to lie in a plane, but this does not appear to be 
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Figure 4. A representation, in the style of Fig. O showing the 
CMB WMAP3 octopole, with the 6*2, (/> distribution with ci = 
0.358 fixed. The third WMAP3 multipole vector with coordinates 
((/), C2) = (48.6° , 0.353) and its antipode are represented as a light 
blue pair of dots. The scale of the plot is the same as in Fig. \3\ 




Figure 5. A further representation of Figs. [3ji and [J] plotted 
now on the surface of the sphere. The multipole vector axes are 
represented by the same colours as before. 



Figure 3. A representation of the 3-dimensional octopole vec- 
tor PDF P3{ci,C2, <f>) from II28I I. in standard configuration. Each 
frame a-j represents a different position of ci < 1 (the green dot 
on the vertical axis), and is a plot over the 62 1 '/'-sphere in Moll- 
weide projection. Frame a has ci = 0, increasing in increments of 
0.1 to frame j, with ci = 0.9 (since P3(1,C2,</') = 0). The fixed 
MV axis in ±2; is represented by red dots at the north and south 
poles, and the opposite vector to vi is shown on each side of the 
plot (at ip = ±7r) as green dots. For each frame (with fixed ci), 
the shading represents the probability P-j{ci, C2, (f>), that is, the 
position of the MV ±tJ2 (parameterized by 02,4')- Contours are 
plotted at values of P3,max/10. If the probability density function 
were uniform over ci , C2 , </>, the uniform density would be about 
the level of the lowest contour line. 



particularly anomalous. The reason for the disagreement 
originates in the particular statistical measure used, and 
various measures have been introduced in the literature, 
based upon their ease for numerical computation. Here, 
we develop an alternative, analytic approach based on 
the full joint probability distribution of the multipole vec- 
tors, which complements the other numerical approaches. 
As can be seen from Fig. |4l the multipole vectors can 
have a significantly planar configuration while not being 
anomalous with respect to a Gaussian distribution. In the 
following subsection, we study octopole planarity further 
by calculating the statistics for a natural measure of MV 
planarity. It should be noted that the discussions of the 
papers cited above, significance is added to the observed 
planarity of the octopole by its correlation with the orien- 



tation of the observed quadrupole; we make no comment 
on inter-multipole orientations here. 

Comparing with the ensemble of the completely inde- 
pendent MVs model, whose probability density is constant 
l/27r ~ 0.159 (about the lowest contour line level in Figs. |31 
[4]) the CMB data is obviously more likely to originate from a 
Gaussian distribution of aim coefficients. The plots of Figs. 
[3^ and 2] are plotted spherically in Fig. [S] 



4.3 Statistics of planarity of ^ = 3 multipole 
vectors 

Using the analytically-calculated probability density func- 
tion for the three multipole vectors of the cosmic octopole, 
it is possible to predict exactly the non-planarity statistics 
of the Gaussian model. 

As we described above, up to a global rotation, the oc- 
topole vectors have three degrees of freedom, and the full 
joint probability density function of these vectors is given 
for the Gaussian model in ((28}. We will compare this with 
the ensemble with completely independent MVs, whose joint 
PDF for ci, C2, is the constant l/27r. 

An analytically natural measure of (non-)planarity of 
the multipole vectors vo,vi,V2, is the magnitude of the 
scalar triple product 



r = \{vo X Vi) ■ W2I 
= |siS2sin0 

for Vi given by (|26|) . where s; 



(30) 

cf. The sign of the 
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Figure 6. Plots of the probability density functions of the mul- 
tipole vector non-planarity r for the two ensembles: Gaussian 
model (black) and uncorrelated multipole vectors (red). 

triple product has no meaning, as the ordering of the vectors 
is arbitrary, as is the choice of ± direction of the Vj. Geo- 
metrically, r is the volume of the parallelepiped with edges 
uo, Hi, V2- Clearly it is if and only if the vectors are copla- 
nar (linearly dependent), and it is 1 if and only if they are 
orthogonal, r does not directly determine the spacing of the 
MVs; in the data, not only are the MVs nearly planar, they 
are also quite uniformly spaced (their angles are all ^ 60°). 
Any planarity measure, including r, depends on the pairwise 
spacing of the MVs as well as their overall planarity. 

In other studies, the normalized angular momentum 
dispersion, maximized over all directions, is used as the 
measure of octop ole planarity (|de Oliveira-Costa et al.|[2004l : 
ICopi et al.ll200 6l). This quantity, which positively weights oc- 
topole configurations with equally spaced MVs (03,3 -1-03,-3 
in some frame), is rather hard to analyze using the analyti- 
cal methods here; r is used as an alternative planarity mea- 
sure which is mathematically natural - (|30|l is the simplest 
quantity to combine three vectors in three dimensions - al- 
though equally spaced planar MVs, as opposed to any other 
planar configuration, are not particularly distinguished, r is 
properly a measure of non-planarity (being smaller the more 
planar the configuration), and we henceforth adopt this de- 
scription. 

In the Appendix, the PDF of r is derived for the Gaus- 
sian model and the independent MVs model. For the Gaus- 
sian model, the PDF is written down as a certain double 
integral which can be integrated numerically for each r. In 
the independent MVs case, the PDF is simply 

Pmd{T) — arccosr. (31) 

Plots of the PDFs for r for the two ensembles are pre- 
sented in Fig. |6l the uncorrelated ensemble tends to have 
a lower (more planar) value of r than the Gaussian model, 
on account of the repulsion of the MVs for the Gaussian 
ensemble. The first two moments in each case are 

{r)i„d = 7v/8 ^ 0.393, 
(r')i„d = 2/9 ^ 0.222, 
('r)Gauss = 0.686, 

{r')GausB = 0.529. (32) 

Thus, for the Gaussian model, although the configuration 
with maximum probability has orthogonal multipole vec- 



0| 




Figure 7. Probability distribution of the multipole vector con- 
figuration conditional on a planar arrangement. The vertical axis 
is 01, the horizontal 82, both in degrees (with ^ Si ^ 90°, 
—90° ^ 62 ^ 90° for all possible planar configurations). 

tors, corresponding to r = 1, the mean r is less than 0.7. 
The data for WMAP3 map, with octopole multipole vectors 
given by (|20p . has a non-planarity measure of 

TWMAP — 0.655, (33) 

which is surprisingly high given the literature on the appar- 
ent anomalous planarity of the octopole; for the Gaussian 
model this value is in very good agreement with the mean. 
As described above, this difference arises in the choice of r, 
rather than other measures, as planarity measure. 

For the third-party maps, r has values 0.749 (TOH), 
0.676 (PARK) and 0.505 (PE). 

The full probability density p8|) may also be used to in- 
vestigate the behaviour of the multipole vectors if they are 
constrained to be planar; in the present parameterization, 
this is equivalent to the restriction (p = (the possibility 
that <j) — TT has been replaced by an extended range for 62, 
—90° ^ 02 ^ 90°). With these restrictions, the probabil- 
ity distribution (for integration with respect to dci dc2) is 
proportional to 

Poau.. (6*1, 1 planar) = (34) 
sin^eisin^6>2sin^(6>i -6>2) 
(9 + cos 26li + cos 26*2 + cos 2(6*1 - 02)y/^ ' 

where A is a normalization constant. This distribution is 
plotted in Fig. [T] Although it is obvious that the most likely 
configuration with enforced planarity is with multipole vec- 
tors maximally distant (i.e. 9i — 60°, 62 = —60°), this dis- 
tribution quantifies this precisely for the Gaussian model. 



5 DISCUSSION 

We have reviewed the multipole vector method of interpret- 
ing CMB maps, including an analytic computation of the 
full probability density function of the multipole vectors for 
the statistically isotropic Gaussian random model. We an- 
alyzed this function for i = 2,3, and a comparison with 
data revealed that the third-year WMAP data is consistent 
with the Gaussian model, despite the apparent anomalous 
planarity of the octopole. 

A natural extension of this work is to provide a full 
PDF analysis of the higher multipoles; however, this is 
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mathematically complicated. P3 depends on three param- 
eters, and therefore can be plotted (as in Fig. |3]) as a 1- 
parameter family of 2-dimensional maps. The next multi- 
pole, £ = 4, requires five parameters to specify the mul- 
tipole vector configuration, and P5 requires 7. In addition 
to this increase of parameters, the term in the denomina- 
tor of the PDF (|15p . (|16|) . summing over {2l)\ permuta- 
tions of the MVs, becomes significantly more difficult to 
compute as I increases. For I > 2Q the problem is exac- 
erbated by the limited numerical accuracy of the various 
MV-finding algorithms. It is clear that as we extend the 
work to higher multipoles, statistics based on derived quan- 
tities, rather than the full PDF, will be necessary. Further- 
more, we ha ve not addressed the importan t issue of inter-^ 
correlations ([Land fc Maguei"ioll2005al . I2OO6I ) , or correlations 
between the CMB and ext er nal re ference frame s - as pointed 
out by Tcgma rk et all (120031): |de O livoi ra-Costa et al.l 
(|2004h : [Schwarz et al.l (|2004h : ICopi eral.1 l|2006l . l2007l ). the 
apparent plane of the octopole is strongly aligned with the 
orientation of the dipole. Such correlations are indicative of 
a breakdown of SI, independent from any assumption about 
the distribution of the a^mS. 

Thus, we defer to a future work a study where measures 
of anisotropy will be constructed from the MVs, in place of 
the full probability density function. These anisotropy mea- 
sures must be defined to be as orthogonal as possible so 
to avoid a repetition of information. Also, the MVs will be 
used to systematically define a Cartesian frame for each mul- 
tipole, so to probe the orientation and inter-.^ correlations of 
the multipoles, and thus probe the issue of SI independently 
of any assumptions about Gaussianity. 

It is fortuitous that the standard model for the CMB, 
namely Gaussian distributed aim coefficients, corresponds 
precisely to t he case studied for spin coh erent SU(2) 
polynomials (Bogo molnv et al] Il996l : lHannavl 1996), with 
the additional constraint ([3]). This means that quantities 
such as the 2-point MV correlation function p2, the full 
MV joint PDF Pt, and, in fa ct, any iV-point MV cor- 
relation function (jDennisI l2005l ) can be calculated exactly 
analytically using the methods of random polynomials. 
This suggests the application to the CMB of other ideas 
from other fields th is technique is used, suc h as random 
spins l|Hannavl Il996l) , Reny i-Werhl entropv |Helling et all 
I2OO6I '). Skyrmion analysis ("Ativah fc Sutcliffel 120021 ) ■ and 
quantum chaos |Bogoniolny ct al. 199^), which may provide 
other powerful tools to analyze cosmological maps. 

One of the most difficult aspects of CMB analysis is, of 
course, that we are provided with only one sample function 
to test our statistics: the observed Universe. This fact is all 
the more apparent at the low multipoles where there are 
only a small number of degrees of freedom. As it is precisely 
these few data points that probe the largest scales in the 
Universe, it is essential that we hone our statistics to exploit 
the low-multipoles of the CMB as carefully as possible. 
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APPENDIX 

Derivation of the full MV PDF (fT5|) 

The coefficients of any polynomial can be written in terms 
of the roots, the so-called symmetric polynomials. In the 
present case. 



(35) 



where s„ is the nth symmetric polynomial for the polyno- 
mial of order 2i (with 2^ + 1 coefficients, n = 0, 21 + 1), 



= n ^ 

l^il<i2<-..<ire^2ffc = l 



(36) 



and the factor aexp ^— jX]i=i 't'i^ is the coefficient for 
Therefore 



Mmflfm = aexp ^— (— 1) 



(37) 



a is a real number since the overall phase factor of ai is fixed. 
Since sq = 1, and S21 = Hiii = exp (2i Yfi=\ 4>i^ , then 



a-e = (— l)^aexp * ^ 0i J S21 = (— l)*a|. (38) 

By hypothesis, the coefficients ao, ai, ae are indepen- 
dent and identically distributed Gaussians (ao is real, the 
others complex). The correctly normalized total probability 
distribution is therefore 



— ;= exp , 



ao 



E 



1 a^js^+fp 



— ;= exp , 

V \ rn — —£ 



1 a' 



^^*+i/2^^Pl^ 2 2(2/)! 

In the second line, |a^mP has been replaced with (|a^m|^ -|- 
|a?.„P)/2 = aWsl+n^\^ + |s|_„p)/2/i^ using JSTj. In the 
final line, the quantity E has been introduced, where 

21 21 
E = ^i!(2£-»)!|s.|^= J2 U^^^ + OCi.)), (40) 

i — cr(zS2ii~^ 

where the second equality, in which E is realized as a 
su m over permut ations a over the 21 roots, was observed 
bv lHannavl (|l996l ). The quantity Ce = YL^-t \aiJ^ I{21 -f 
1), described in Section E] equals a^E/2(2/)! by ((SS]). The 
PDF of any completely random ensemble depends only on 
this quantity. 

In order to write the probability density (|39p as a dis- 
tribution over the multipole vector root configuration, it is 
necessary to find the appropriate Jacobian {a^m} ^ {Ci}- 



(39) 
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Using standard methods (e.g. iBogomolnv et al.l l|l996l )) it is 

straightforward to show that 



V 



T\m=-i Ia*™] n,=i iCil 



-dad^^{G}. (41) 



where D is the square root modulus of the polynomial dis- 
criminant, 



V ■ 



n 1^'-^^- 



(42) 



i.j — l.i<j 



The final joint PDF for the multipole vectors is the in- 
tegral of this quantity with respect to a (—00 < a < 00) 
(using the Gaussian probability density (I39|l . It must be di- 
vided through by (.\ (since the £ independent roots are in- 
distinguishable), and 2^ (since Ci)Ci+< ^re indistinguishable 
for j = l,...,^): 



V 



da a exp 



(£-1/2)! 



2(2£)!' 

(2£)K+l/2-p 



(43) 



Rearranging the numerical prefactor, {£ — l/2)!/7r^+^''^ = 
(2^-l)!/(27^)^ and {2£y.'+'/y UL=-,M = U-U jl Fi- 
nally, substituting E via (|40p gives (|16|l . 

In order to transform from the stereographic plane to 
the direction sphere, it is necessary to multiply by the stere- 
ographic Jacobian factor 2~^f|*^j(l -I- ICiP)'^ (this factor 
differs from the usual stereographic Jacobian factor by 2^, 
since the labelling MVs are in the northern hemisphere, 
meaning that ^ 1 for j = 1, ...,£). Substituting the 
roots by (|14l) . with appropriate interpretation on the direc- 
tion sphere f Han nav.1996, ) gives the joint probability density 
function (1151). 



Derivation of the r PDF 

The probability density function for r is found by substitu- 
tion of (/) — ^ r by 



(44) 



The Jacobian of this transformation is 
d(l) 1 



dr 



(45) 



The derivation of the PDF for r in the independent MVs 
model is easy, as the PDF is constant with respect to the 
volume element dci dc2 dep. The PDF can therefore be calcu- 
lated from the ci,C2,(j) normalization integral, with </!> sub- 
stituted by r using (|44)| and (|45l) . then integrating ci and 
C2. This integral is 



2 /■! /■! {1-4) a -4) 

1 — — I dci I dc I dr 



V(l-c2)(l-c?)- 



(46) 

Without changing the integrand, the order of integration 
is changed to C2, then ci, then r last, with ^ C2 ^ 



\/l - rV(l - cl), ^ ci s; 1 - then sC r sC 1. Fi- 
nally, ci and C2 are rescaled in order that the volume of 
integration is again a cube: 



1 pi rl 



1 — — dr dci / dc2 



Jo Jo 



(47) 

ci and C2 may now be integrated directly. 



1 = / dr arccosT, 
'0 



(48) 



giving the PDF of equation (|3ip . 

The strategy for the Gaussian model is identical, al- 
though the integrand in (I46p is multiplied by the Gaussian 
model PDF ((29}, with substituted by r. The final Gaus- 
sian r PDF therefore can be written in terms of a double 
integral. 



„ , , 1500^37? , f\ 

-PGauss(T) = / dCl I dS2X 



7!" JO J-1 

4/2 2 2/ , \2\/ 2 2 2 2\ 

a (a —T] Si[ciC2ri + as2) )[a — r/ S1C2) 



(49) 



C2 {r/^sKScl + {ci_C2T) + as2)(ciC2r? -|- 3as2)) + a2(5 -I- Zclr)'^)) 
where 



r) = Vl si = \/l - c^, C2 = si, a = 1 - vi'^cl. 

(50) 

As in (|34p and Figure [T] negative values of S2 appear as they 
replace values of (f) between 7r/2 and 37r/2. 
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